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Abstract: We show that the conventional threshold resummation calculation cannot de- 
scribe well the low energy Drell-Yan (DY) data without including the non-perturbative 
correction terms which are deduced from analyzing the asymptotic behavior of the resum- 
mation formalism. It is demonstrated that the non-perturbative correction is generally 
small for the large invariant mass DY pairs produced at the Tevatron and the LHC. 
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1. Introduction 

With ample data accumulated at the Fermilab Tevatron and expected at the up-coming 
CERN Large Hadron Collider (LHC), precision measurements at hadron colliders become 
possible. To test the Standard Model (SM) predictions and to probe possible new physics 
signatures, many higher order theoretical calculations have been performed in recent years 
to study the phenomenology of Drell-Yan (DY) pair, W boson, top quark pair and Higgs 
boson produced in hadron collisions. It is well known that in perturbative QCD (pQCD) 
hadronic cross sections could receive logarithmically enhanced contribution in the partonic 
threshold region, which corresponds to z ^ 1, where z = Q'^/s with Q being the invariant 
mass of the final state (say, the DY pair) and \/I the center-of-mass (CM) energy of 
the partonic process. These contributions take the form of "plus" distributions, such as 
[ln"'(l — z)/{l — z)]^, and can be resummed to all orders in the expansion of the strong 
coupling as ■ This is called the threshold resummation calculation which plays an important 
role for precision tests at hadron colliders. 

A typical Sudakov exponent in threshold resummation in the modified minimal 
subtraction (MS) ^ |3| scheme, 

Ef^{N)=- [\z'-^^:^^\ C ^A(/»)[a,(AQ2)]+^(/0[«^((l_,)2g2)]l 
Jo I- Z ^ J{l-z)2 A J 

generally involves singularities, known as the Landau pole of QCD running coupling. This 
occurs when the relevant energy scale is smaller than the QCD scale Aqcd, as z — > 1 
or A — > 0, where the pQCD fails and the non-perturbative QCD effects must set in. To 
avoid directly confronting the non-perturbative contribution, a few approaches have been 
proposed in the literature. One approach is to apply the principal value resummation Q 
to choose a contour to evaluate the above integral without hitting Landau pole |5|. 
Another method to treat the Landau pole problem was recently proposed by Bonvini et al. 
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in Ref. p. Other approaches were also proposed in the hterature. For example, in Ref. Q, 
Becher et al. established an approach to resum the logarithmic contributions directly in 
the momentum space based on effective theory. In Ref. Kidonakis et al. proposed to 
approximate the resummation calculation by an expansion of the resummed cross section, 
which will be further discussed below. In this paper, we propose a new approach to 
add non-perturbative correction terms to the minimal prescription threshold resummation 
formalism |^], as suggested by the perturbative expansion of joint resummation. 

This paper is organized as follows. In Section 2 we shall first briefly review the relevant 
part of threshold resummation (RES) formalism. In Section 3 from the threshold resum- 
mation formalism we induce the functional form of the non-perturbative (NP) correction 
terms to be added to the minimal prescription threshold resummation calculation. Then, in 
Section 4 we apply this new threshold resummation formalism, denoted here as RES+NP, 
to three DY experiments with CM energy of the hadron colliders ranging from 38.76 GeV 
to 1960 GeV [^, 10, 11 1. Prom the comparison to data, we determine the coefficients in 
the NP correction terms of the threshold resummation formalism. We also compare the 
RES+NP result with the predictions of several popular approaches. Finally, we give our 
prediction for the DY pair production at the LHC. Section 5 contains a brief summary of 
the conclusions. 



2. Threshold Resummation 

The differential cross section of the DY pair, in terms of its invariant mass (Q) and rapidity 
(y), in the hadron CM frame (with energy \/5) is I! 



dQ^dY Jq 7o ^1 Jo X2 xiX2 dzdy 

(2.1) 

where r = Q'^/S. y is the rapidity of lepton pair in the parton CM frame, f^"- are the 
parton distribution functions (PDFs), and fip the factorization scale. 

Conventionally, the rapidity-integrated cross sections only need the Mellin transform to 
turn the convolution into multiplication. However, for the case of the rapidity distribution 
the Fourier transform with respect to Y is also needed [|^, |l^ . Applying the Mellin- Fourier 
transform to Eq.(|2.1|), we obtain 



rl f-ln^ ^ H1H2 

Jo JinJT dQMY 



dQ^dY^ (2.2) 

pi \- 



where 



and 



f^^iN ± iM/2) = [' dxx^-i±^^^/2/pT'(^0 (2-3) 
Jo 

P'P^{N,M)= dzz^-^ dye'^'y%— (2.4) 

Jo Jlny^ dzdy 
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are the PDFs and differential cross section in moment space, respectively. As shown in Ref. 
1 13], near threshold region the M dependence in liJ^'^'^^{N, M) is negligible, and it is a good 
approximation to take the well-known resummed form |l^ for the rapidity-integrated 
cross section. 



5?i^^(iV,M) ^ ^?^r{N) = exp \Y^E^^'{N)] exp [ ^^2 / ^^^^/^,^(a,(/.2)) 



xiv|iJPiP2(a,(4))Pexp 
/Of 



Q/N 



X SP'PHl,a 



'■\N2 



Pexp 



(2.5) 



TP1P2 

^ s 



where N = N exp (je), 1e is the Euler constant, and /z/? is the renormalization scale. P 
denotes path ordering The first exponent in Eq.(^]^), as given in Eq. (|1.1|), resums 
the collinear and soft gluon emission contributions from the initial state, with 



Aif^\as) = C,(^ + \K(^ 



vr 



TT 



(- 


--) 


V18 


6 ) 




vr 



(2.6) 

(2.7) 
(2.8) 



where nj is the flavor number of light quarks. Cf = Cp = {N^ — l)/{2Nc) for initial 
state quarks and Cf = Ca = -^c for initial state gluons. Nc is the number of colors. In 
Eq.(^^), at the one-loop order, = (as/7r)(|Ci? — Cplri-N) for quarks, and jg/g = 
(Qs/7r)(/3o — C^lniV) for gluons. The /? function is defined as 



1 dings 
2^ di^ 



n=0 



a,\("+2) 



TT 



with 



(5o = {IICa - 2ny)/12, 
(5i = {17C\ - BCAUf - 3CFnf)/24. 



(2.9) 

(2.10) 
(2.11) 



Ts in Eq. (|2.q ) is the soft anomalous dimension matrix, which can be derived from the 
eikonal diagrams j^, |3| and is given at the one-loop order by 



Ts = (1 - vrz). 



TT 



At the next-to-leading logarithm (NLL) accuracy, w^^p^ approximated as [13| 
u^^'lliN) = I^jP^'P'Cia,) exp[5i(A) IniV + g^iX)], 

with 

A = /Jqcks In A^/tt, 



(2.12) 

(2.13) 
(2.14) 
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2A + (1 - 2A)ln(l - 2A) 



(2.15) 



52(A) 



2A + ln(l - 2A) + - ln2(l - 2A) 



+ ^[2A + ln(l-2A)]ln 



^[2A + ln(l-2A)] 



(2.16) 



2Aln 



2 ' 



where lv^^^ is the Born differential cross section in moment space. By matching the 



moments of the NLO cross section, the coefficient C{as) in Eg. ( 2. 13 ) can be obtained as 



follows fl^ 



TT V 3 



2tt'' 3, g2 
- In —75- 
2 



(2.17) 



The cross section in momentum space can be obtained via inverse Mellin-Fourier transform 



1 /"C+^oo -| foo 

'P\z,y) = / dNz~^— dMe-'^-^yuP'P^{N,M). (2.18) 

27rz Jc-ioo 27r J^^ 



In order to avoid double-counting the fixed order contributions up to the next-to-leading 
order (NLO), it is needed to subtract the first two orders of Os expansion to obtain the 
final cross section, i.e.. 



HiH2 
RES 



,da 



H1H2 

NLO 



+ S 



da 



H1H2 

NLL 



Hi Ha 
NLL 



dQ^dY dQ^dY dQ^dY \ dQ^dY J as=o \das dQ^dY J as=o 



a., 



a da^^^^ 

.5 NLL 



, (2.19) 



where o-Sf^^, corresponding to cj^tV^q in the momentum space, denotes the NLO cross 



section. 



For comparison, in this paper we also investigate the next-to-next-to-leading order 
(NNLO) expansion |8|, |l^] of the resummed cross sections. At first, in order to reproduce 
the fixed order expressions, we recover the M dependence in 



^PiP2 M) = - + + (AT - iM/2)] , (2.20) 



although we have neglected the M dependence in the numerical evaluation of the resummed 
cross section [13|. For the inverse Mellin-Fourier transform of the logarithms ln^{N±iM/2) 
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we find 

1 fC+ico 1 

' dNz-""— / dMe-'^y 



27ri Jc-ioo 27r 
=5(1 - z)6{y) 

--5{l-z)6{y±Unz) 

1 i-C+ioo -I 

--— dNz~''5{y± -In z), (2.21) 

2v^^ Jc-ioo 2 

1 fC+ico 1 



/ diVz-^ — / dMe-'^y [In^ {N ± iM/2) + 7^5] 

Jc-ioo 27r 
1 fC+ioo 1 /-oo 

= ^ / dN'z-""'— dMe-^^^(^±5i--)[ln'=(Ar')+^^] 
27rz yc-ioo 27r y^^^ 

1 f-C+ioo 1 

= / diVz-^ln'=(iV)(5(y ±-lnz). (2.22) 

2Tn Jc-ioo 2 

Hence, we can replace the inverse Fourier transform and M dependence by the function 
[5(y + ilnz) + 5(y - ilnz)]/2. The remaining calculation is the same as the rapidity- 
integrated differential cross section. If we expand the exponent of the resummed cross 
section in Eq.(^) to the NLO accuracy []l6[ , the differential cross section for DY assuming 
HF = fJ'R = Q 



where Dk{z) = [ln'=(l - z)/{l - z)]+, C3 = 4Cf and ci = 2CfC2 - ^Cf- wg'^' is the 
coefficient of 5(1 — z)5{y) in the Born differential cross section uj^'^^ . It is evident that this 
reproduces the dominant contribution near threshold in the NLO differential cross section, 
as given in Ref. Similarly, at the NNLO accuracy, the expansion of the exponent in 

Eq. ( p^) yields 



^^e^^'^'\z,y) =a;r'(^)'{ic2Z)3(z) - PoCsD^iz) + (C3C1 - + 2CfK)D^{z) 

+ (Csc^ - 2(3^ci)Do{z) + (^c? + icici - \ucl)5{l - z)] (2.24) 
5{y+\\D.z)+5{y - \ \ryz) 

X . 

2 

Below, we define the NNLO-NLL (next-to-next-to-leading order and next-to-leading loga- 
rithmic) [^8| corrected differential cross section as 

, ,PlP2 _ , ,PlP2 I , ,PiP2(2) (n 9cr\ 

'^NNLO-NLL — ^NLO ^ '^exp. ■ t^Z.ZOj 

3. Nonperturbative Effect 

While deriving Eq. (|2.13| ) from Eq.( |l.lD , approximation ^ has been made to avoid Landau 
pole in the original expansion, and the integration in Eq.(|l.l|) was carried out via per- 
turbative expansion. This approximation could fail if the NP contribution in Eq.(l.l) is 
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large. Hence, we propose to add NP correction terms in the resummation formalism to 
better approximate the total contribution from Eq.(^]^). To find out the proper functional 
form to parameterize the NP corrections, we examine the joint resummation formalism in 
Ref. where the resummed cross section in moment (and impact parameter) space is 
expressed as 

(3.1) 



Jeik) 



(eik) , 



where 



2NkT 



^ ^ i=a,b 



+ In 



(3.2) 



Z)(f)(iV,6) 



dkiA, 



+ 



abyoisikr 



, kj 



In 



kx 

'q 



+ e 



-ib-kn 



Kn 



2NkT\ 
Q ) 



dk 



T 



dk'^Wab{k'^,kT,as{np)) 







(3.3) 



-ib-kx 



[Ko(2N^ 



Ikl + k^ 



Ko( 



2Nh 



T 



V Q 



+ In 



'4 + fc2 



The infrared renormalon singularities occur as A;t ^ for a large Q value with Q ^ Nkx- 
Below, we shall ignore the 6-dependent term which is only relevant to transverse momentum 
resummation. From the expansion of the Bessel function K(){x) for small value of argument 
X, K(){x) ~ — ln(xe'''^/2), we find that the leading behavior of E^^^^^ in fey ^ limit can 
be described by the following two functions: 



iV2 , iV2 

and —r In , 
Q2 Q2 vivgo 



(3.4) 



Since Aab behaves as 1/Q^ when /c-r — > with a large Q value ||T^, d'^^^^^ behaves as 
which is suppressed as compared to E^^^^\ Another correction term at the order of 



but suppressed by l/N'^, can be obtained by examining the behavior of e'^j^^^ and D 
in the limit that /ct — > and Q is not very large as compared to /ct, i.e., Q ~ Nkx- Since 
for large argument x, Kq{x) ~ e~^/y^, we can neglect the contribution from Kq. The 



^(eik) 
ab 



leading contribution in this limit comes from the logarithm term In^kx/Q) in D^^' 
implies that 

J-lnf^ 

should be considered. 

After combining the above three sources, we obtain the NP correction term 



[eik) 



which 



(3.5) 



Ar2 



ai + 02 In 



Q 



NQo 



1 



+ 03— 2 In 



Q_ 

Qo 



(3.6) 
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where the parameter Qq is chosen to be 1.3 GeV^ and the dimensional NP parameters 
ai, 02 and 03 are to be determined by comparing the theory prediction of RES+NP with 
experimental data. In this improved threshold resummation formalism, cD^^p^ Eq.(^ 
is written in the form 



^'re!lMN) = «^(A^)exp[5^P(Ar)] 



(3.7) 



4. Numerical Results 

In our numerical calculation, the SM parameters are chosen as in Ref. 



The running 



QCD coupling is evaluated at the three-loop order |2C], and CTEQ6.6M PDFs |21] are used 
with = fJ-R = Q. In Table ||, we list the set of DY experimental data to be considered in 
our analysis. We follow the method of least chi-square (x^) analysis in Ref. |22] to find the 



Experiment 


no. of data points 


\/5(GeV) 




E605 g 


119 


38.76 


15% 


E866 (pp) 


184 


38.76 


6.5% 


cDFig 


29 


1960 


5% 



Table 1: The DY data sets considered in the analysis, cr^ is the experimental normalization 
uncertainty of the nth experiment. 

best fit by allowing the overall normalization of each experiment to vary. They are denoted 
as Ni, N2 and Y3 for E605, E866 and CDF Run-2 DY (via Z/j* production) experiments, 
respectively. The NP parameters are determined by the best fit of the RES+NP calculation 
to the set of DY data, which yields 



01 = -0.60 GeV^ 



02 



-2.87 GeV^ 03 = 5.58 GeV^ 



with Ni being around 1 and P^r degree of freedom (dof) 1.07, cf. Table ^. 

In Table ^ , we also compare the result of RES+NP calculation with a few other theory 
calculations which include the NLO [|l2|, 0, |2|, NNLO |2|], NNLO-NLL |§, and the usual 
threshold resummation calculation (RES) |jl^. . We repeat the same fitting procedure for 
each theory calculation by allowing the normalization of each data set to float within its 
experimental uncertainty in order to find the best fit to theory prediction. 





X|605(^l) 


X|866(^2) 


X^Df(^^3) 


Xtotal 


xVdof 


NLO 


103 (0.99) 


247 (0.97) 


19.5 (0.95) 


369.5 


1.11 


NNLO 


271 (1.25) 


772 (1.28) 


20.6 (0.92) 


1063.6 


3.20 


NNLO-NLL 


147 (1.04) 


908 (1.05) 


17.4 (0.94) 


1072.4 


3.23 


RES 


196 (1.18) 


1198 (1.20) 


17.6 (0.97) 


1411.6 


4.25 


RES+NP 


128 (1.03) 


209 (0.96) 


17.0 (0.97) 


354.0 


1.07 



Table 2: The minimal value of and normalization factor for each experiment and theory pre- 
diction. 



""^This is the energy scale that the CTEQ6.6 parton distribution functions are evolved from. 



As shown in the table, the NLO results agree well with the data because the E605 data 
were included in determining the CTEQ6.6M PDFs at the NLO. The NNLO result is about 
a factor 3 worse than the NLO result in x^/dof, owing to the low energy E605 and E866 
data, which indicates that a NNLO PDF set is needed to improve the comparison. The 
NNLO-NLL prediction is similar to the RES results. The conventional threshold resumma- 
tion calculation (RES) cannot describe well the E605 and E866 data with the mass of the 
DY pair ranging from 7 GeV to 18 GeV and 4.2 GeV to around 15 GeV, respectively. The 
largest difference between the results of RES and RES+NP occurs in the E866 data. To 
examine it in more detail, we show in Fig.|^ the comparison among various theory calcula- 
tions for one particular set of E866 data, with 0.60 < xp < 0.65, as an example, where xp 



is the Feynman-x variable [10|. To clearly examine the difference between the theoretical 
predictions and the experimental data, we define x = {Dn — Tn/N2)/a^, where Dn, cr^ and 
Tn denote the data value, experimental measurement uncertainty and the theoretical value 
for the nth data point, respectively. As shown in the figure, the largest deviation from data 
occurs when Q is around 5 to 6 GeV. The RES result also fails to describe data unless the 
NP correction terms are included which is the result of RES+NP. Hence, we conclude that 
to describe the low energy DY data with the threshold resummation formalism, the NP 
correction terms must be included in order to take into account the part of contribution 



missing from approximating the Sudakov integral Eq.(Ll) by its perturbative expansion 
Eq.(Pl. 



^ • . 0.60<Xp<0.65 


= ^ . . • * ^ ^ X 

— X X X « « 

~ .... I .... \ .... \ .... I .. . 


NNLO-NLL 
X NNLO 

* RES 

• RES+NP 

1 .... 1 .... 1 


« Si 

1 



> 

O 



X 
T3 
O 
T3 



10' 

















• E866 






NNLO 






- - - RES 

— RES+NP 






. . 1 .... 1 





9 10 

Q (GeV) 



11 



12 



13 



14 



Figure 1: Different theoretical predictions compared to the E866 (pp) experimental data, with 
0.60 < xf < 0.65. In the lower plot, the theoretical results have been multiplied by I/N2. 

Next, let us compare various theory calculations to the large Q (around 100 GeV) DY 
data, taken by the CDF Collaboration at the Tevatron Run-2. 

As shown in Fig. |2|, for Y > 2, the NLO result becomes smaller than the data while the 
NNLO, RES and RES+NP calculations give similar results, though the RES+NP result 
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Figure 2: The ratio of CDF Run-2 DY data to various theoretical predictions (after being muhi- 
phed by I/N3). The shaded area indicates the statistical error of the data. 



gives the lowest this set of data, cf. Table |2[ We could also compare the total 

cross section of the DY pair produced at the Tevatron Run-2 and the LHC. The result of 
comparison is listed in Table ^. It shows that the results of RES and RES+NP are about 



a{pp/pp l^l + X)(pb) 


NLO 


NNLO 


RES 


RES+NP 


Tevatron 


240.7 


242.0 


248.0 


247.98 


LHC 


2047.9 


2036.8 


2115.2 


2115.3 



Table 3: The total cross sections of the DY pair (with 66GeV < Q < 116GeV) produced at the 
Tevatron Run-2 and the LHC. 



the same, and differ from NNLO (and NLO) by about 2.5% at the Tevatron and 3.7% at the 
LHC. Hence, we conclude that the effect from the NP correction terms to the conventional 
threshold resummation formalism for large Q values in high energy hadron collisions is not 
important. Finally, we show in Fig. ^ the differential cross section as a function of the 
invariant mass of the DY pair produced at the LHC. In this figure, dasvu/dQ represents 
the differential cross section including the RES-I-NP result and the electroweak box diagram 
contribution, originated from the re-scattering of W and Z bosons in loop diagrams |25]. In 
order to study in detail the higher order correction to the shape of the da/dQ distribution, 
we also plot the ratio of differential cross sections normalized by dasvu/dQ. The label 
RES-I-NP indicates the ratio of da jies+np /dasuu, etc. It is evident that the results of 
RES and RES-I-NP are almost the same except when the value of Q is small, less than 
about 30 GeV. 
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5. Conclusion 



In conclusion, we have investigated the NP effect in the threshold resummation calcula- 
tions. From analyzing the asymptotic behavior of the conventional threshold resummation 
formalism (RES), we proposed the NP correction terms to be included in the improved 
threshold resummation formalism (RES+NP). They are subsequently determined by fit- 
ting to the DY data from E605, E866 (pp) and CDF Run-2 experiments. We found that to 
describe the low energy DY data with relatively small value of invariant mass, the NP effect 
cannot be ignored in threshold resummation calculation. In contrast, the minimal prescrip- 
tion threshold resummation formalism (RES) gives a similar prediction as RES+NP when 
the invariant mass of the DY pair is large. Though it remains to be seen how well the 
threshold resummation formalism could describe the top quark pair production rates at 
the Tevatron and the LHC |2C], the type of NP effect discussed in this paper is not expected 
to be important. On the contrary, this effect will become relevant for describing the low 
invariant mass DY pairs produced at the Relativistic Heavy Ion Collider (RHIC) at the 
Brookhaven Laboratory. 
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